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SUMMARY 


The  conducted  research  was  directed  at  broadening  the  science  base  for  LEC 
type  melt  growth  of  InP  through  the  establishment  of  quantifiable  thermal  boundary 
conditions  in  the  melt  and  through  the  stabilization  of  the  melt  by  means  of  applied 
axial  magnetic  fields.  In  the  extension,  this  research  was  aimed  at  systems 
development  which  have  direct  applicability  in  crystal  growth  technology. 

In  view  of  the  non-availability  of  high  temperature  heat  pipes  throughout  the 
duration  of  this  research  program,  much  of  the  work  was  carried  out  with  analog 
model  systems.  This  work  yielded  semiconductor  single  crystals  in  which  all 
micro-segregation  inhomogeneities  -  which  characterize  melt  grown  semiconductors  - 
as  well  as  radial  macro-segregation  effects  were  completely  suppressed.  Related 
modelling  efforts  revealed  that  pressure  effects  substantially  influence  the  growth  and 
segregation  behavior  as  well  as  the  generation  and  distribution  density  of  dislocations 
during  conventional,  non-encapsulated  growth.  Pressures  up  to  35  atmospheres  were 
not  found  to  significantly  influence  segregation  during  LEC  growth  with  B203  at  layer 
thicknesses  of  1  cm  and  more. 

The  fundamental  concepts  of  this  research  effort  were  the  basis  for  the  most 
advanced  research  program  on  heat  and  mass  transfer  control  with  magnetic  melt 
stabilization  during  low  pressure  LEC  growth  of  GaAs  sponsored  jointly  by  DARPA/Air 
Force  Materials  Laboratory  and  for  the  extended  research  phase  on  "Development  of 
Model  Based  Magnetic  LP-LEC  Growth:  Large  Diameter  GaAs.”  Published  results 
obtained  in  conjunction  with  this  research  have  stimulated  industrial  development 
efforts  in  the  area  of  high  temperature  heat  pipes.  Most  recently 
Wah-Chang/Teledyne  has  announced  the  development  of  a  prototype  device  which 
will  be  tested  at  MIT  and  applied  to  growth  of  GaAs  and  InP. 
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I.  Segregation  in  RF  Heated  High  Pressure  LEC  growth 


Continuing  in  an  effort  to  elucidate  the  growth  behavior  during  high  pressure 
LEC  growth  of  InP,  a  study  of  the  segregation  behavior  was  undertaken.  This  study's 
aim  was  to  clarify  differences,  if  any,  between  high  pressure,  RF  heated  Czochralski 
growth  and  conventional  low  pressure  resistance  heated  growth  of  semiconductors  in 
moderately  dimensioned  crucibles  (6  cm  diameter)  with  and  without  liquid 
encapsulation.  The  model  system  Ge:Ga  was  chosen  for  the  following  reasons: 

(1)  quantitative  segregation  data  is  easily  obtained  from  the  elemental 
system  Ge:Ga; 

(2)  the  melting  point  similarity  between  Ge  and  InP  makes  it  an 
appropriate  model  for  general  segregation  characteristics  expected 
in  the  growth  of  InP; 

(3)  the  equilibrium  segregation  coefficient  of  Ga  in  Ge  is  comparable  to 
Sn  (a  common  dopant  used  in  producing  n  type  material)  in  InP. 

This  study  was  undertaken  to  investigate  the  difference,  if  any,  between: 

(1)  RF  and  resistance  heating, 

(2)  encapsulated  vs.  unencapsuiated  growth,  as  well  as 

(3)  low  vs.  high  ambient  pressure. 

It  was  expected  that  axial  micro-segregation  in  RF  heated  systems  would  be 
be  much  more  severe  than  that  observed  in  resistance  heated  growth  due  to  the 
unavoidable  problem  of  allignment  of  the  susceptor  within  the  RF  coil  and  the 
differential  coupling  that  would  be  experienced.  This  effect  was  expected  to  lead  to 
pronounced  thermal  asymmetry  in  the  growth  enviroment. 

Experimental  Results 

A  series  of  Ge:Ga  crystals  were  grown  from  quartz  crucibles  under  conditions 
in  which  most  of  the  critical  growth  parameters  were  systematically  varied. 

Subsequent  to  growth,  crystal  sections  containing  the  rotational  axis  of  the  growing 
crystal  were  cut  with  a  diamond  saw.  These  samples  were  subsequently  (1)  mounted 
on  brass  disks  with  ln:Sn  solder  to  form  back  side  ohmic  contacts.  (2)  iapped  with 
5  ixm  Al203,  and  (3)  polished  with  a  solution  of  200  ml  Nalcoag  2355  colliodal  silica, 
400  ml  Dl  water,  and  6  mi  H202. 

After  achieving  a  mirror  finish  on  the  sample  surface,  the  crystals  were  etched 
in  a  solution  of  HF:CH3C00H:H202  (1:1:1)  to  reveal  compositional  variations  (growth 
striations).  These  samples  were  then  subjected  to  single  point  spreading  resistance 
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measurements  to  obtain  quantitative  data  related  to  the  segregation  behavior  of  Ga  in 
the  Ge.  Focus  was  placed  upon: 

(1)  axial  macro-segregation:  the  determination  of  the  effective 
segregation  coefficient; 

(2)  axial  micro-segregation: 

(3)  radial  macro-segregation  to  obtain  information  related  to  the  shape  of 
the  solute  boundary  layer,  as  well  as  information  related  to  bulk  melt 
flows. 

The  matrix  of  growth  experiments  undertaken  is  listed  in  table  1. 

Results  and  Discussion 

The  Effect  of  Pressure  and  Encapsulation  on  Segregation 

Crystals  were  grown  to  investigate  the  effects  of  pressure  and  encapsulation 
upon  segregation  in  the  standard  high  pressure  growth  enviroment.  Typical 
micro-segregation  data  from  these  crystals  is  presented  in  fig.  1.  (Note  that  the 
non-encapsulated  crystal  turned  polycrystalline  immediately  upon  a  pressure  change; 
hence,  there  is  no  data  for  this  portion  of  the  crystal.)  The  charge  carrier  profile  is 
essentially  flat,  showing  a  lack  of  any  evidence  of  rotational  striations.  Previously,  it 
was  reported  that  encapsulation  in  small  scale  systems  had  the  effect  of  eliminating 
micro-segregational  inhomogeneities  by  creating  a  transition  from  turbulent  melt 
convection  to  laminar  forced  convection  through  the  establishment  of  (a)  a  reduced 
destabilizing  thermal  gradient  and  (b)  retarding  the  convective  melt  flows  through  the 
establishment  of  a  shear  boundary  layer  at  the  melt/encapsulant  interface.  In  the 
present  study  on  growth  with  larger  crucible  sizes,  it  is  seen  that  the  encapsulant  does 
not  bring  about  these  transitions,  as  indicated  by  the  turbulence  dominated 
micro-segregation  profile.  This  can  be  understood  through  consideration  of  the 
increased  Reynolds  number  with  increasing  dimensions,  and  hence  a  larger  driving 
force  for  turbulent  thermal  natural  convection.  The  radial  segregation  data  exhibit  a 
flat  profile,  indicative  of  either  a  highly  uniform  or  virtually  non-existent  solute  boundary 
layer,  consistent  with  a  highly  turbulent  melt. 

The  axial  macro-segregation  data  allowed  for  the  calculation  of  k  erf  based 
upon  the  Scheil  equation: 
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From  axial  macro-segregation  data,  the  effective  segregation  coefficient  for 
the  high  and  low  pressure  regions  of  sample  1  were  calculated.  These  values  were 
found  to  be: 


4 


*58 


low  pressure  (50  psi): 


ken  =  062 


high  pressure  (100  psi):  k«tt  =  0  60 


Considering  the  error  limits  in  the  measurement  technique,  it  can  be  stated  that 
the  two  values  are  the  same;  hence  the  change  in  pressure  had  no  apparent  effect 
on  the  macroscopic  segregation  behavior  It  is  noted  though,  that  these  values  are 
extremely  high  for  the  Ge.Ga  system  (k0  =  0.087).  An  analysis  of  the  boundary  layer 
thickness  for  this  experimental  finding  (k«»  )  according  to  the  BPS  theory  yields: 
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R  =  growth  rate  =  1.11  x  1 0  "3  cm/s 
D  =  1 .9  x  10_4  cm2  /s 
kQ  =  0.087 
Cs  /  C,  -  0.6 

5  =0.41  cm 


When  comparing  this  8  value,  complications  in  the  data  interpretation  become 
evident  with  that  obtained  through  a  Cochran  analysis  for  a  rotating  disk; 

1/3  1/6  -1/2 

8  =  1 .6  D  v  oj 

where: 

D  =  1.9  x  102  cm  /s 
v  =  1 .4  x  1 03  cm  ^  /s 
o)  =  rotation  rate 
8  =5.2x1 0-^  cm 


It  is  observed  that  the  BPS  analysis  gives  a  boundary  layer  thickness  which  is  in 
total  contradiction  with  the  establishment  of  convective  melt  flows.  It  is  also  noted 
that  for  much  of  the  growth,  the  boundary  layer  thickness  as  predicted  by  BPS  is 
found  to  be  greater  than  depth  of  the  melt. 
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Resistance  vs.  RF  Heated  Growth 

To  clarify  the  anomalously  high  k  err  as  well  as  the  reasons  for  the  lack  of 
rotational  striations  (periodic  compositonal  variations),  several  crystals  were  grown. 

The  axial  micro-segregation  data  from  a  crystal  grown  in  an  LP-Resistance  heated 
furnace  (fig.  2)  shows  periodic  compositonal  fluctuation  which  correlated  with  the 
seed  rotation  rate.  There  is  no  evidence  of  non-rotational  compositional  variations. 
This  lack  of  turbulent  convection  in  the  case  of  resistance  heated  growth  can  be 
explained  by  considering  the  fundamental  differences  between  the  heat  transfer  in  the 
two  modes  of  heating  (fig.  3).  In  the  case  of  RF  growth  (with  a  conventional 
susceptor  configuration)  a  large  driving  force  for  convection  is  expected  due  to  the 
general  bottom  heating  nature  of  the  thermal  configuration. 

Radial  etch  and  charge  carrier  profiles  (fig.  4)  are  seen  to  exhibit  an  essentially 
flat  growth  morphology  consistent  with  a  uniform  solute  boundary  layer.  Axial 
macro-segregation  data  again  suggest  a  k  etr  of  approximately  0.6  (figs.  5  &  6). 

Anomalous  Axial  Macro-Segregation 

In  response  to  the  anomalously  high  k  etf  found  under  all  growth  conditions,  a 
series  of  crystals  were  examined  to  clearify  this  effect.  Related  representative  axial 
macro-segregation  data  is  presented  in  figures  7-11.  In  fig.  7  the  Scheil 
equation's  graphical  analysis  is  shown.  It  is  seen  that  this  analysis  leads  to  a  keff  of 
0.519.  In  figures  8  and  9,  the  concentration  data  is  plotted  along  with  the  theoretical 
concentration  profiles  for  k  eti  equal  to  0.1  and  0.519  respectively.  These  graphs 
show  that  the  concentration  data  follow  a  k  eft  of  -  0.1  up  to  g  =  0.5,  while  following 
a  keff  of  -0.5  for  the  second  half.  Figures  10  and  11  show  raw  data;  in  this  case, 
the  effective  segregation  coefficient  is  in  the  range  expected  for  conventional  Cz 
growth.  This  inconsistent  behavior  indicates  that  a  condition  of  uncontrolled  bulk  flow 
is  affecting  the  segregation  behavior  in  a  manner  which  makes  reproducible  growth 
and  hence  comparison  of  data  impossible. 

Conclusions 

1.  The  heat  transfer  differences  between  conventional  RF  heating  and  resistance 
heating  gives  rise  to  a  large  driving  force  for  turbulent  convection  in  RF  heated 
growth. 

2.  Encapsulation  does  not  bring  about  a  transition  from  turbulent  to  laminar  flow  in 
a  growth  system  of  moderate  dimensions. 

3.  Pressure  has  no  detectable  effect  upon  the  macro  and  micro  segregation 
characteristics. 
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Sample  Furnace  Type 


1 

2 

3 


Table  1 


Encapsulant  SRR 


HP-RF 

No 

HP-RF 

Yes  (1 

cm) 

HP-RF 

Yes  (1 

cm) 

LP-Res. 

No 

LP-Res 

NO 

PR 


4  cm/hr 
4 


Pressure 


50/100  psi 
50/100  psi 

300/500  psi 
1  atm 
1  atm 

100  psi 
1  atm 
1  atm 
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II.  Heat  Transfer  In  High  Pressure  Growth  of  InP 


Analytical  calculations  of  the  relative  contributions  of  radiational  and  convective 
heat  transfer  from  the  surface  of  an  InP  crystal  grown  at  ambient  pressures  from 
0-40  atm.  are  presented.  These  calculations  were  then  used  as  a  guide  to  the 
modification  of  a  finite  element  heat  transfer  model  of  LEC  GaAs  growth  for  InP 
growth. 

Excess  thermal  stress  has  been  considered  to  be  one  of  the  major  causes  for 
dislocation  generation  in  III— V  compound  semiconductors  [1-3].  In  thermal  stress/heat 
transfer  models  of  LEC  growth,  focus  is  on  GaAs.  These  analyses  are  often  adapted 
to  InP,  but  the  effect  of  the  high  pressure  ambient  gas  is  neglected  in  heat  transfer 
calculations.  For  this  reason,  the  convective  contribution  to  heat  transfer  relative  to 
radiational  heat  transfer  has  been  analytically  calculated.  These  calculations  were  then 
used  as  a  guide  in  the  modification  of  the  model  of  Motakef  [3]  to  take  into 
consideration  the  physical  properties  of  InP  as  well  as  the  inclusion  of  geometric  and 
convective  heat  transfer  considerations  specific  to  the  growth  of  InP. 

Radiation  vs.  Convection:  The  Effect  of  the  Ambient  Pressure 

To  evaluate  the  relative  contributions  of  radiation  and  convection  to  heat 
transfer  as  a  function  of  pressure,  the  non-dimensionalized  heat  flux  from  the  crystal 
surface  for  each  effect  was  calculated.  The  radiational  flux  is  formulated  as: 

q-  =  ea  F(T?-T*) 

R 

q  FT  R  n(  ®  1  “  ®2  ) 


where: 
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radiation  number 
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The  non-dimensional  convective  heat  flux  is  expressed  as: 

q'c=  h  (T,  -  l  ) 

~q;=  Bi  (0t  -  ©2  ) 


where: 


Bi  = 


where  : 


h  =  0.525  Ra 


“{i} 


=  0.525  (Pr  Gr)  1/4(k/D) 

Using  the  above  equations,  the  relative  contributions  to  heat  transfer  were 
calculated  for  ambient  pressures  between  1  and  40  atmospheres  and  temperature 
differences  (between  the  crystal  surface  and  the  ambient)  of  400  and  600  degrees. 
Table  1  shows  the  calculated  contributions  which  are  presented  graphically  in  fig.  1  as 
( q  conv  /q  rad )  vs.  pressure.  It  is  seen  that  at  the  pressures  used  for  the  LEC  growth 
of  InP  (P  >  27.5  atm),  convective  heat  transfer  contributes  more  than  1/4  to  the 
possible  heat  transfer  in  the  system.  This  graph  suggests  that  the  effects  of 
convective  heat  transfer  are  negligible  in  modelling  the  growth  of  low  pressure 
systems,  but  convective  heat  transfer  must  be  considered  in  high  pressure  systems. 


Finite  Element  Modelling 

Constant  crystal  diameter,  absence  of  crystal  shoulder 

As  the  first  attempt  in  modelling,  the  finite  element  code  of  Motakef  [3]  was 
modified  to  account  for  the  thermo-physical  properties  of  InP  (table  2)  and  to  allow 
for  convective  heat  transfer  boundary  conditions.  In  this  model  a  cylinderical  crystal 
was  assumed  to  grow  from  a  melt  encapsulated  with  an  optically  transparent 
encapsulant.  Two  cases  were  studied:  one  without  and  one  with  convective  heat 
transfer  boundary  conditions  at  the  crystal/gas  interface  as  well  as  at  the 
encapsulant/gas  interface.  Figures  2  and  3  show  the  excess  thermal  stress  profiles 
at  the  crystal  center,  half  radius  and  the  periphery  as  a  function  of  non-dimensional 
axial  location.  From  these  two  plots  it  is  clear  that  convection  in  the  gas  phase  alters 
both  the  positions  and  the  magnitude  of  the  maximum  excess  stress.  The  stress 
level  in  the  crystal  was  found  to  increase  by  about  25%  in  a  high  pressure 
environment.  This  indicates  that  any  heat  transfer  control  approach  in  the  growth  of 
InP  must  address  convection  in  the  gas  phase. 
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Crystal  with  conventional  shoulder  and  a  transparent  or  opaque  encapsulant 


Since  a  cylinder  does  not  closely  model  the  actual  shape  of  Czochralski  grown 
InP,  it  was  decided  to  evaluate  the  effect  of  shoulder  geometry  upon  the  temperature 
distribution  in  the  crystal.  In  this  case  it  is  important  to  consider  that  non-linearities  in 
the  temperature  field  give  rise  to  thermal  stress.  Figure  4  presents  the  calculated 
temperature  field  for  an  InP  crystal  with  a  60°  shoulder  encapsulated  by  a  transparent 
encapsulant.  The  isotherm  profile  exhibits  pronounced  curvature  in  the  shoulder 
region;  this  finding  can  be  interpreted  as: 

1.  The  shoulder  acts  as  a  funnel  through  which  the  heat  must  move  creating 
a  restriction  to  heat  flow  and  hence  extreme  curvature  in  the  isotherm 
profile; 

2.  The  shoulder  surface  radiates  heat  to  the  cold  upper  furnace  wall,  thus 
drastically  increasing  the  heat  transfer  from  this  surface. 

In  the  case  of  an  opaque  encapsulant,  an  inflection  in  the  isotherm  curvature 
can  be  noted  (fig.  5).  The  magnitude  of  curvature  is  less  than  for  the  transparent 
encapsulant  which  suggests  that  a  thermally  opaque  encapsulant  may  be  more 
favorable  for  thermal  stress  reduction  than  the  semi-transparent  B203  used 
extensively.  It  must,  however,  also  be  considered  that  a  thermally  opaque 
encapsulant  will,  through  its  isolating  effect,  substantially  reduce  radial  thermal 
gradients  in  the  melt  and  thus  impede  diameter  control  during  growth. 

Shoulder  and  opaque  encapsulant  at  two  stages  during  growth 

Since  the  dislocation  generation  and  density  distribution  will  be  affected  by  all 
of  the  stress  distributions  that  the  crystal  experiences  during  growth,  two  different 
stages  of  growth  were  modeled,  both  while  the  crystal  was  still  being  greatly  affected 
by  the  shoulder  region.  In  the  first  case  (fig.  6)  the  crystal  is  fully  immersed  in  an 
opaque  encapsulant.  As  seen  in  fig.  7,  the  stress  distribution,  already  at  this  early 
stage  of  growth,  greatly  exceeds  CRSS.  The  stress  distributions  at  varying  axial 
locations  (figs.  8-12)  show  the  existence  of  the  now  classical  “W"  shaped  radial 
profile.  At  most  locations  away  from  the  crystal  melt  interface  the  excess  stress 
profile  is  found  to  be  relatively  low  and  flat. 

A  later  stage  in  growth  (after  the  shoulder  has  emerged  from  the  encapsulant 
but  the  crystal  has  not  been  rounded  to  constant  diameter)  is  depicted  in  fig.  13.  In 
this  case,  we  see  that  the  crystal  experiences  two  maxima  in  its  stress  profile  (fig. 
14),  one  at  the  crystal/melt  interface  and  one  at  the  encapsulant/gas  interface. 
Examining  the  radial  profiles  of  the  stress  distribution,  a  much  more  pronounced  “W" 
shape  in  the  stress  distribution  as  compared  to  earlier  stages  of  growth  is  observed. 
As  with  the  earlier  case  (figs.  15-21),  the  level  of  stress  decreases  and  the  degree 
of  variation  decreases  in  the  axial  direction  above  the  encapsulant.  This  is  in 
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qualitative  agreement  with  the  dislocation  density  distribution  (both  radial  and  axial) 

observed  on  InP  crystals  grown  in  this  laboratory. 

Conclusions 

1.  The  high  pressure  ambient  gas  has  a  significant  effect  on  the  heat  loss  from  the 
crystal  surface,  contributing  as  much  as  35%  to  the  overall  heat  transfer. 

2.  The  geometry  of  the  crystal  (the  shoulder)  serves  to  increase  the  lack  of 
planarity  in  the  temperature  profile  in  the  crystal,  giving  rise  to  larger  thermal 
stresses. 

3.  Even  during  the  earliest  stages  of  growth,  when  the  crystal  shoulder  is  still  fully  in 
the  encapsulant,  the  matrix  experiences  severe  excess  stress,  suggesting  that  in 
the  case  of  InP  fully  encapsulated  growth  likely  may  not  be  effective  in  reducing 
the  dislocation  density  in  InP. 
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Thermoelastic  and  Growth  Properties 

of  InP 

InP 

Thermal  Conductivity  (W/cm  K) 

1215  T"1  324 

Suface  emissivity 

0.5 

CRSS  (  dynes/cm  K) 

•1  q(1  24  +  58.05/T) 

Young's  Modulus 

9.51  x  1011 

Poison's  Ratio 

0.287 

Elasticity 

7.5  x  10  ~6 

Encapsulant  (Boric  Oxide) 

Thermal  Conductivity  (W/cm  K) 

2.37  x  10-3  +  1.1  : 

Growth  System 

Ambient  Gas:  Argon 

Temperature: 

700  K 

Viscosity 

5.72  x  10  ~  4 

Thermal  conductivity 

0.044 

Furnace  Wall  Temperature: 

400  K 

Encapsulant  Radiation  Number 

3.29 

Crystal  Radiation  Number 

0.0758 

Trans  i.  1  (With  Convection) 
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III.  NIR  Microscopy  of  N-type  III— V  Compound  Semiconductors 

Compositional  inhomogeneities  have  a  major  effect  upon  the  electrical  and  op¬ 
tical  properties  of  semiconductor  materials.  Techniques  which  permit  the  qualitative 
and  quantitative  characterization  of  these  variations  (on  a  microscale)  in  III— V  semi¬ 
conductors  are  virtually  non-existent.  Near  infrared  microscopy  provides  a  technique 
which  shows  the  potential  of  filling  this  need. 

Free  Carrier  Absorption 

N-type  GaAs  and  InP  exhibit  free  carrier  absorption  in  the  infrared.  This  ab¬ 
sorption  phenomenon  holds  a  potential  for  compositional  characterization  of  these  ma¬ 
terials.  Free  carrier  absorption  is  a  non-direct  process  involving  the  scattering  of 
electrons  in  the  conduction  band.  In  examining  the  potential  of  this  effect  for  materi¬ 
als  characterization,  the  theoretical  calculations  of  Walukiewicz  et  al.  [1,2]  are  used 
as  a  basis. 

In  modelling  free  carrier  absorption  in  GaAs  and  InP  Walukiewicz  considered  the 
effect  of  acoustical  phonon  scattering,  optical  phonon  scattering,  and  ionized  impurity 
scattering  upon  the  absorption  coefficient.  The  total  absorption  due  to  free  carriers 
thus  becomes: 


at(  X)  =  aop(  \)  +  o^c(  \  )  +  Ojmpf  M 

where: 


a  op,\1 

a  op.  \0  I 

f  x  1  ' 
^  X  0  ' 

)2'5 
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0 

(  K  i  N 

<  X  Q  / 

r 

(  ^  1  ' 

.  35 

aimp,x1  = 

aimp,  Xq  1 

X  x  0  ' 

) 

In  figures  1  and  3  the  theoretical  free  carrier  absorption  characteristics  of 
GaAs  and  InP  with  electron  concentrations  of  1016/cm3  and  10,8/cm3  are  presented. 
From  these  curves  it  can  be  seen  that  the  absorption  coefficient  is  largest  in  the 
mid-IR,  with  the  absorption  coefficient  of  InP  rising  more  rapidly  than  that  of  GaAs.  In 
figures  2  and  4,  the  same  data  is  plotted  in  the  wavelength  range  of  interest  to  near 


infrared  microscopy.  It  can  be  seen  from  these  two  graphs  that  a  moderately  doped 
sample  will  show  very  little  wavelength  dependence  of  the  absorption  coefficient.  In 
fact,  this  doping  range  would  require  thick  samples  for  any  measurable  transmittance 
change  to  occur.  The  transmittance  is  given  by: 

T  _  (1  ~  R)2  exp(-CLd) 

1  -  R^exp(-2ad) 

where: 

R  =  Reflectance 
T  =  Transmittance 
d  =  sample  thickness 
a  =  absorption  coefficient 

The  general  validity  of  these  calculations  can  be  seen  in  fig.  5  where  the  free 
carrier  absorption  of  lnP:Sn  crystal  is  shown. 

For  near  IR  transmission  light  microscopy  to  be  a  viable  approach  to  micro¬ 
segregation  characterization,  the  absorption  coefficient  must  show  sensitivity  to  the 
free  electron  concentration.  Figure  6  shows  this  dependence  for  a  radiation  with  a 
wavelength  of  1  pm  In  excellent  agreement  with  previously  reported  absorption 
measurements  on  n-type  GaAs  [3],  a  nearly  linear  relationship  between  the  absorp¬ 
tion  coefficient  and  the  electron  concentration  is  exhibited. 

Experimental  Conditions  for  Analysis 

NIR  characterization  requires  an  optical  system  consisting  of  a  uniform,  stable 
light  source  and  a  sensitive  detecting  device.  Also  required  is  the  ability  to  quantify 
the  transmitted  light  intensity  sensed  by  the  detector.  To  accomplish  this,  an  Axitron 
microscope  (Zeiss)  in  combination  with  a  silicon  based,  extended  range  infrared  TV 
camera  (Dage-MTI)  were  used.  The  resulting  image  was  processed  utilizing  a 
512  x  512  x  256  image  processing  system  (Recognition  Technology)  resident  in  a 
super-mini-computer  (MassComp).  This  system  allows  for  the  quantification  of  the 
light  intensity  at  every  pixel  position  across  the  field  of  view.  Thus,  theoretically  the 
absorption  coefficient  and  electron  (dopant)  concentration  can  be  determined  for  any 
position  in  a  given  wafer. 

The  potential  of  this  approach  is  shown  in  fig.  7,  in  which  a  near  infrared 
micrograph  on  a  GaAs: Si  sample  is  shown.  With  calibration  this  distribution  will  be 
translated  into  a  real-time  measurement  of  the  dopant  distribution. 

Fourier  Analysis 

Utilizing  a  light  intensity  distribution,  information  related  to  the  spatial  distribution 
of  inhomogeneities  can  be  obtained.  Through  Fourier  analysis  of  this  distribution,  the 


harmonic  relationships  among  the  various  features  were  calculated.  In  this  analysis, 
the  light  intensity  is  expressed  as: 

...  1  v  i  co  x 

l(x)  =  *2  2,  c  n  e  n 

where: 

1  r  -  icon  x 

cn  =  —  J  l(x)  e  dx 


L 

where  L  is  the  fundamental  period  of  the  function  and  Cn(con)  is  the  "Fourier  spec¬ 
trum"  of  the  function  l(x).  In  fig.  8  a  near  infrared  micrograph  of  an  lnP:Sn  crystal  is 
shown  with  a  super-imposed  light  intensity  distribution  and  Fourier  spectrum.  In  this 
material  the  effects  of  turbulent  melt  convection  are  reflected  in  the  random  nature  of 
the  observed  striae  as  well  as  in  the  resulting  Fourier  spectrum  which  shows  numerous 
modes  indicating  the  non-periodic  nature  of  the  dopant  distribution. 
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